On tumor development: fractional transport approach 
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A growth of malignant neoplasm is considered as a fractional transport approach. We suggested 
that the main process of the tumor development through a lymphatic net is fractional transport of 
cells. In the framework of this fractional kinetics we were able to show that the mean size of main 
growth is due to subdiffusion, while the appearance of metaphases is determined by superdiffusion. 

PACS numbers: 87.17.-d, 05.40.Fb 



I. INTRODUCTION 



Application of fractional calculus in biology is mainly 
used to describe an activity of living systeni(see recent 
reviews [l|). Our purpose is to describe a tumor growth 
process by means of cell transport inside a fractional net- 
work. In the present consideration we propose a link 
between tumor spread and fractional transport, whose 
mathematical apparatus is well established Q, y, |j{ We 
propose simplified mathematical models using heuristic 
arguments on tumor growth. On the basis of statement, 
we are citing '^: "Lymphatic spread is more typical of 
carcinomas whereas the hematogenous route is favored 
by sarcomas. There are numerous interconnections, how- 
ever, between the lymphatic and vascular systems, and 
so all forms of cancer may disseminate through either 
or both systems." - we focus on a possible mechanism 
of this growth as a result of fractional transport of can- 
cer cells along lymphatic net system, which has a fractal 
Hausdorff dimension [l| . Using a simplified approach of 
fractional transport, we are also suggesting a possible an- 
swer on a question how neoplasm cells appear arbitrary 
far from a main tumor. The tumor growth is a compli- 
cated process which consists of consequences of various 
and different phases. This phases can be considered as 
independent processes with different time duration. The 
last depends on various environmental conditions as well. 

The growth of mahgnant neoplasm can be described 
as a two-step process P- The first phase is the growth 
and differentiation of the primary tumor and the second 
phase can be described as a metastatic phase, during 
which there are many changes in the primary tumor's 
cells that lead to the metastasis. The difference between 
the methastatic cells and the cells of the primary tumor 
includes, for example, a different expression of adhesion 
molecules such as E-cadherins, catenins etc. j^]. Those 
differences are discussed in many studies, see e.g. Q, 
according to which after certain proliferation of the pri- 
mary tumor the metastatic stage begins - the cells start 
to express different molecules. Those cells separate from 
the primary tumor - a process is named metastasis. 



The metastasis cells can be spread throughout the lym- 
phatic net system or throughout the blood stream (vas- 
cular system) ^. Apparently, there are cytokines which 
support the lymphatic proliferation and encourage the di- 
rectional movement toward lymphatic system |8|. There 
are tumors which express different specific chemokine re- 
ceptors and they are more likely to spread through the 
lymphatic system. 

In the present paper we suggest a simplified scheme of 
tumor spread along the lymphatic net. In our scheme we 
are introducing two - step process, where the first one is 
a biological process of cells fission. The second phase is 
cells transport itself inside lymphatic net structure. We 
primary focus on the second phase, supposing that the 
main size of the tumor development is due to fractional 
transport, and present the two times tumor development 
model (TTT-model). 



II. MAIN ASSUMPTIONS FOR THE TWO 
TIMES TUMOR DEVELOPMENT MODEL 

In order to construct the TTT-model for the tumor 
development description we are introducing a number of 
assumptions which are essential for our approach: 

(a) We suppose, that there are two time scales Tf and 
%. The first one corresponds to the volume grow due 
to proliferation of cells by fission. During this time cells 
are strongly interacted [3- There is no cell transport, 
and only the fed cells can proliferate. The duration of 
Tf could be arbitrary large and it can reach 10^ — lO^sec 

5]. The second time % corresponds to cell transport. 
We suppose that during this time there are enough cells 
without fission, and the interaction between tumor cells 
is rather small. Hence the cells form an initial packet of 
free particles, which spread along the lymphatic net. 

(b) The lymphatic net is a system of organic tubes with 
complicated crosses and arbitrary (random) number of 
entrances at nodes. It fulfilled by lymphatic liquid with 
laminar flow. Citing ji5j,|10j "About 100 milliliters per our 
of lymph flows through the thoratic duct of a resting hu- 
man and approximately another 20 milliliters flows into 



the circulation each hour through other channels, making 
a total estimated lymph flow of about 12Qml/hr, that is 
between 2 and 3 liters per day" . The Hausdorff dimension 
of this geometry is fractal d^ < 3 (3 is embedding dimen- 
sion) . Since the geometrical dimension of the lymphatic 
net is fractional, the neoplasm development corresponds 
to anomalous diffusion. 

(c) The tumor spreading process consists of the following 
time consequences 



Tfil)T,{2)Tfi3). 



(1) 



There are different realizations of this chain of times, due 
to different duration oiTf(i) and Tt{i), where i = 1, 2, . . . 
and it means simply a step-number of the process, ft 
should be underlined that the processes with different 
steps i are absolutely independent and there is possible 
any realization of the time duration Tt{i). Therefore, 
one can introduce the probability distribution function 
(PDF) for Tt{i)- As a result of these realizations, there 
are few scenarios of the tumor development. This pro- 
cess depends on the rate of the tumor cell spreading. 
If transport is subdiffusive, that is typical for the lym- 
phatic spread of carcinomas j^], it can be described in 
the framework of continuous time random walk (CTRW) 
0, [l3| ■ The transport can be superdiffusive due to so- 
called Levy-type flights when cells may traverse all of 
the lymph nodes ultimately to reach the vascular com- 
partment via the thoratic duct. In this case for any 
anatomic localizations metastasis are possible. It should 
be stressed that an every cell (a particle) carries with it- 
self its own trap, and this is the principal deference from 
the standard CTRW, where traps are external with re- 
spect to the transporting particles. 



III. SUBDIFFUSIVE GROWTH 

There are many possibilities to describe the cell trans- 
port as a diffusive process taking place during the time 
scales Tt{i) in the framework of the fractional calculus. 
One of the simplest approximation for this process is 
CTRW. We used here a simplified model which neglects 
proliferation of cells. It means that a number of cells 
participated in the transport is conserved. Suppose that 
Tf{i) > "^(j), and following the CTRW construction 
(see, for example, Q), we believe that during %{i) a 
cell "jumps" with the jump length variance 



a^ = 



l^^{l,t)dldt 



(2) 



and Tf{i) are waiting times elapsing between two succes- 
sive jumps, where 



r 



t-ip{l,t)dldt 



(3) 



the waiting time distribution has the following form 

/•OC 
m^ / i^{l,t)dlr.l/{l+t/Tf + ", (4) 

Jo 

with < a < 1 the behavior for a CTRW is subdiffusive 
with the mean squared displacement being 



{f{t))^e. 



(5) 



is the mean waiting time. These values are characterized 
by a probability distribution function (PDF) ij}{l,t). If 



This subdiffusive relaxation process is responsible for the 
size of the primary tumor. It may be relevant if one 
takes into account that during times Tt{i) there is diffu- 
sion on fractal lymphatic net of dimension df . The mean 
squared displacement in this case is given by (jSJ, where 
the transport exponent a depends now on df. In these 
cases the relaxation process is modeled by fractional dif- 
fusion equation, where the PDF ip[l, t) has non-Gaussian 
form. The square root from mean squared displacement 
(P) is an average size of the main tumor development. 

When a > 1, the transport on the fractal lymphatic 
net is superdiffusive. It may describe a metastasis pro- 
cess. In this case \/{l'^{t)) is an average distance between 
a new tumor location and the primary tumor after time 
t. This time i ~ T is a time when a cell moves long 
enough without fission. 



IV. COMB MODEL 

In what follows we consider a simple model of frac- 
tional transport of neoplasm cells. In the framework of 
this model we are able to determine the size of the tu- 
mor development. Anomalous diffusion on a comb model 
is an oversimplified presentation of diffusion on the frac- 
tal lymphatic net. The model was proposed in [l3J for 
drift particles in percolation systems and used for some 
applications in solid state jlj| and biophysics [l3- The 
comb-like geometry shown in Fig. 1 is a toy model of a 
porous diffusion medium. We borrow this consideration 
from [iJI primary to describe subdiffusion of cancer cells 
on the lymphatic net and demonstrate a fractal nature 
of cancer cells transport as well as to determine a size of 
cancer growth. A special feature of diffusion in this model 
is a possibility of motion in the x-direction only along the 
structure axis, when y ^ 0. In this case motion in the 
y-direction is a trap, where a particle stays in the comb 
teeth for some time before it moves along the structure 
a;-axis. This model approximates cell transport along the 
lymphatic net embedded in 3D sphere {R, 6, (/>). Diffusion 
along the structure axis (y=0) specifies motion along the 
lymphatic net leading to the increasing of the radial size 
R = R{t) of the tumor. Motion along the y-direction 
corresponds to cell transport in the polar and azimuthal 
directions. This motion does not lead to the increase of 
the radial size of the tumor and may be considered as a 
trap. The PDF of waiting or delay times for the infinite 
teeth is ip{t) ^ t~^ , where 1 < /? < 2, such that the mean 
waiting time is diverged. It follows from the presented 



FIG. 1: The comb structure 



geometry, that the diffusion coefficient D^x differs from 
zero only when y = 0, namely D^x = Di5{y). The dif- 
fusion coefficient along the y direction is Dyy = D2 and 
Dxy = Dyx = 0. Therefore, The Fokker-Planck (FP 
equation for the Green function G = G{x, y, t) reads 



FP) 

M 



f-*)0 



D. 



9y2 



- 6{x)6{y)6{t) . 



(6) 



In the Laplace-Fourier (p, k) space after carrying out cor- 
responding transformations for the x, t variables we ob- 
tain from Eq. (jBJ the following solution 14] 



G{k,y,p) 



e^^[-{p/D2f/^\y\ 
2{D2pY'^ + Dik^ 



(7) 



For diffusion along the structure axis one can infer the 
Green function c(a;,t) — G{x,0,t) 14]. Therefore, the 
mean squared displacement of mobile particles is anoma- 
lous subdiffusion: 



{x^it)) 



{x^{t)c{x,t)) 
{c{x,t)) 



= Dl{TTt/D2 



,1/2 



(8) 



This mean squared displacement {x'^it)) is a mean radial 
size of the tumor after elapsing transport time t. For the 
TTT model it means that a general elapsed time is 



N 



t = tN^Y.^r}{i) + %{i)) = tf + u. 



(9) 



During the fission time tf neoplasm cells do not move, 
then t — tt \n ^ . If tt is large enough, this process may 
give a possible explanation not only for the mean size of 
the primary tumor but also it could be mean distance 
between the primary tumor and a metastasis. 

The transport (or grow) process accelerated by includ- 
ing so-called Levy-type flights into consideration. In the 



framework of the comb model it still leads to subdiffusion 
with transport exponent 1/2 < a < 1 [T3. Therefore to 
obtain a superdiffusion process one needs to go beyond 
the comb model consideration. The transport can be 
superdiffusive due to Levy-type flights when cells may 
traverse all of the lymph nodes ultimately to reach the 
vascular compartment via the thoratic duct. In this case 
the Levy-type process is a result of the cell transport in 
the vascular system. Therefore, the solution of Eq. ^ 
should be modified in the following way: along the struc- 
ture axis we add Levy flights with a Levy distribution for 
the jump length with the asymptotic behavior |j] 



ip{x) ~ |a; 



-i-p 



(10) 



where 1 < /i < 2. The solution Eq. {T)) is modified, for 
example, in the following form \\J\ 



GLik,0,p) 



1 



D2P^/^ + Di\k\t^ ' 
The power-law asymptotics for the Green function is 



(11) 



c{x,t) 



|i+M 



(12) 



with g^ > 0. For example, when fi — I, g^j^ ^ 1. In this 
case the mean squared displacement diverges: 



{x\t)) 



(13) 



This expression of Eq. (|13|l could be a simple explanation 
of appearance of metaphases at any long distance from 
the primary tumor. 



V. CONCLUSION 

A cancer growth process is considered in the framework 
of the TTT model. It is shown that the main process of 
the tumor development is a fractional transport of cancer 
cells on the lymphatic net. In the framework of this frac- 
tional kinetics we were able to show that the mean size 
of the main growth is due to subdiffusion, while the ap- 
pearance of metaphases is determined by superdiffusion. 
The underlying mechanism of the tumor development is a 
partition of this grow on fission times Tf(i) and transport 
times Tt{i). Due to this assumption, there are free cancer 
cells, which contributes to fractional transport. There- 
fore, both the characteristic size of the cancer growth 
and metaphases are related to the fractional transport 
exponents a, /i, which depends on geometrical fractional 
dimension of the lymphatic net df . 
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